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Abstract: This work presents results of ab-initio simulations of continuous 
wave transport in disordered absorbing waveguides. Wave interference ef- 
fects cause deviations from diffusive picture of wave transport and make the 
diffusion coefficient position- and absorption-dependent. As a consequence, 
the true limit of a zero diffusion coefficient is never reached in an absorbing 
random medium of infinite size, instead, the diffusion coefficient saturates 
at some finite constant value. Transition to this absorption-limited diffusion 
exhibits a universality which can be captured within the framework of the 
self-consistent theory (SCT) of localization. The results of this work (i) 
justify use of SCT in analyses of experiments in localized regime, provided 
that absorption is not weak; (ii) open the possibility of diffusive description 
of wave transport in the saturation regime even when localization effects 
are strong. 

© 2013 Optical Society of America 

OCIS codes: (290.4210) Multiple scattering; Coherence and statistical optics: (030.1670) Co- 
herent optical effects; Materials: (160.2710) Inhomogeneous optical media. 



References and links 

1. P. M. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw-Hill, New York, 1953). 

2. M. C. van Rossum and T. M. Nieuwenhuizen, "Multiple scattering of classical waves: microscopy, mesoscopy, 
and diffusion," Rev. Mod. Phys. 71, 313-371 (1999). 

3. S. Chandresekhar, Radiative transfer (Dover, New York, 1960). 

4. L. V. Wang and H. Wu, Biomedical Optics: Principles and Imaging (Wiley-Interscience, 2007). 

5. D. Vollhardt and P. Wolfle, "Diagrammatic, self-consistent treatment of the anderson localization problem in 
d < 2 dimensions," Phys. Rev. B 22, 4666-4679 (1980). 

6. J. Kroha, C. M. Soukoulis, and P. Wolfle, "Localization of classical waves in a random medium: A self-consistent 
theory," Phys. Rev. B 47, 11093-1 1096 (1993). 

7. P. W. Anderson, "Absence of diffusion in certain random lattices," Phys. Rev. 109, 1492-1505 (1958). 

8. B. A. van Tiggelen, A. Lagendijk, and D. S. Wiersma, "Reflection and transmission of waves near the localization 
threshold," Phys. Rev. Lett. 84, 4333^336 (2000). 

9. N. Cherroret and S. E. Skipetrov, "Microscopic derivation of self-consistent equations of anderson localization 
in a disordered medium of finite size," Phys. Rev. E 77, 046608 (2008). 

10. C. Tian, "Supersymmetric field theory of local light diffusion in semi-infinite media," Phys. Rev. B 77, 064205 
(2008). 

11. C. Tian, arxiv p. 1301.6225. 

12. M. Storzer, P. Gross, C. Aegerter, and G. Maret, "Observation of the critical regime near anderson localization 
of light," Phys. Rev. Lett. 96, 063904 (2006). 

13. H. Hu, A. Strybulevych, J. H. Page, S. E. Skipetrov, and B. A. van Tiggelen, "Localization of ultrasound in a 
three-dimensional elastic network," Nat. Phys. 4, 945-948 (2008). 

14. Z. Q. Zhang, A. A. Chabanov, S. K. Cheung, C. H. Wong, and A. Z. Genack, "Dynamics of localized waves: 
Pulsed microwave transmissions in quasi-one-dimensional media," Phys. Rev. B 79, 144203 (2009). 

15. T. Sperling, W. Buhrer, C. M. Aegerter, and G. Maret, "Direct determination of the transition to localization of 
light in three dimensions," Nat. Phot. 7, 48-52 (2013). 



16. B. Payne, A. Yamilov, and S. E. Skipetrov, "Anderson localization as position-dependent diffusion in disordered 
waveguides," Phys. Rev. B 82, 024205 (2010). 

17. B. Payne, T. Mahler, and A. Yamilov, "Effect of evanescent channels on position-dependent diffusion in disor- 
dered waveguides," Waves in Random and Complex Media 23, 43-55 (2013). 

18. C. Tian, S. Cheung, and Z. Zhang, "Local diffusion theory for localized waves in open media," Phys. Rev. Lett. 
105, 263905 (2010). 

19. J. B. Pendry, "Symmetry and transport of waves in one-dimensional disordered systems," Adv. Phys. 43, 461-542 
(1994). 

20. A. A. Chabanov, M. Stoytchev, and A. Z. Genack, "Statistical signatures of photon localization," Nature 404, 
850-853 (2000). 

21. J. Wang and A. Z. Genack, "Transport through modes in random media," Nature 471, 345-348 (2011). 

22. E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, "Scaling theory of localization: 
Absence of quantum diffusion in two dimensions," Phys. Rev. Lett. 42, 673-676 (1979). 

23. L. I. Deych, A. Yamilov, and A. A. Lisyansky, "Scaling in one-dimensional localized absorbing systems," 
Phys. Rev. B 64, 024201 (2001). 

24. L. Sapienza, H. Thyrrestrup, S. Stobbe, P. D. Garcia, S. Smolka, and P. Lodahl, "Cavity quantum electrodynamics 
with anderson-localized modes," Science 327, 1352-1355 (2010). 

25. C. W. Beenakker, "Random-matrix theory of quantum transport," Rev. Mod. Phys. 69, 731-808 (1997). 

26. P. W. Brouwer, "Transmission through a many-channel random waveguide with absorption," Phys. Rev. B 57, 
10526-10536 (1998). 

27. A. Yamilov, R. Sarma, B. Redding, B. Payne, H. Noh, and H. Cao, "Position-dependent diffusion of light in 
disordered waveguides," unpublished (2013). 

28. A. Mirlin, "Statistics of energy levels and eigen-functions in disordered systems," Phys. Rep. 326, 259-382 
(2000). 

29. E. Akkermans and G. Montambaux, Mesoscopic Physics of Electrons and Photons (Cambridge University Press, 
Cambridge, UK, 2007). 

30. S. John, "Electromagnetic absorption in a disordered medium near a photon mobility edge," Phys. Rev. Lett. 53, 
2169-2172 (1984). 

31. A. Lagendijk, B. van Tiggelen, and D. S. Wiersma, "Fifty years of anderson localization," Phys. Today 62, 24-29 
(2009). 

32. B. Payne, J. Andreasen, H. Cao, and A. Yamilov, "Relation between transmission and energy stored in random 
media with gain," Phys. Rev. B 82, 104204 (2010). 

33. P. A. Lee, D. A. Stone, and H. Fukuyamak, "Universal conductance fluctuations in metals: Effects of finite 
temperature, interactions, and magnetic field," Phys. Rev. B 35, 1039-1070 (1987). 

34. P. W. Anderson, D. J. Thouless, E. Abrahams, and D. S. Fisher, "New method for a scaling theory of localization," 
Phys. Rev. B 22, 3519-3526 (1980). 

35. B. L. Altshuler, P. A. Lee, and R. A. Webb, eds., Mesoscopic Phenomena in Solids (North Holland, Amsterdam, 
1991). 

36. V. D. Freilikher, M. Pustilnik, and I. Yurkevich, "Effect of absorption on the wave transport in the strong local- 
ization regime," Phys. Rev. Lett. 73, 810-813 (1994). 

37. A. Yamilov and B. Payne, "Classification of regimes of wave transport in quasi-one-dimensional non- 
conservative random media," J. Mod. Opt. 57, 1916-1921 (2010). 

38. L. Zhao and C. Tian, unpublished (2013). 



1. Introduction 

The diffusive description of electromagnetic wave propagation in inhomogeneous media JUG] 
is overwhelmingly successful with wide ranging applications from astrophysics [3 1 to biomed- 
ical optics J4J. Wave interference effects are usually negligible because the phases of the mul- 
tiply scattered partial waves are random. However, for a wave that returns to its source location 
there is always the time-reversed path which necessarily yields the same phase. When the return 
probability is sufficiently high, the constructive interference between the wave propagated via 
the two paths can suppress long-range transport so that the diffusion coefficient in an infinitely 
large (passive) system turns to zero [5 6 1. This is Anderson localization - the phenomenon first 
conceived in context of electronic (de Broglie) wave transport in condensed matter physics Q. 

The question of how the diffusion coefficient evolves with an increase of system size from 
the unrenomalized value of Do toward the limit of D = has been addressed in Ref. flUED 
by extending the original self-consistent theory of Vollhardt and Wolfie |5]|6| to systems of 



finite size. The key prediction of the modified self-consistent theory (SCT) is that the diffusion 
coefficient is no longer a constant but varies spatially. This conclusion is also reached inde- 
pendently in the supersymmetric field theory iflO) . The position dependence of the diffusion 
coefficient arises because the return probability responsible for the renormalization of diffu- 
sion is position-dependent in a system of finite size. As the wave explores the larger and larger 
neighborhood of a source point, the return probability becomes sensitive to the proximity of a 
boundary where the wave has a chance to escape from the system. In fact, Ref. 0T| points out 
that waves propagation in open media can be thought as a highly unconventional macroscopic 
diffusive phenomenon. 

Although SCT has been successful in interpreting the results of several experiments (|T2j- 
[151 . the key prediction of position-dependent diffusion has not been verified in an experiment 
because it requires the access to the inside of the random medium. Position-dependent diffusion 
has been also predicted by the super-symmetry approach [10], suggesting that this is more 
than just a convenient mathematical abstraction. Indeed, quantitative agreement between SCT 
and numerical simulation has been found in disordered waveguides at the onset of localization 
II16II171 . Numerical simulations in Ref. fTSl , show that deeper in the localization regime, the 
prediction by SCT for the position-dependent diffusion coefficient is quantitatively incorrect. 
Instead, by using the supersymmetric field theory it was found that the diffusion coefficient 
exhibits novel scaling lfl8l . This is because SCT underestimates the energy density inside the 
bulk of medium that is strongly affected by the presence of necklace states |fT9l formed via 
resonant tunneling. 

Material absorption and other sources of dissipation are usually a nuisance in experiments 
with classical waves because it leads to the same exponential scaling of e.g. optical conduc- 
tance g exp(— L/^ a o) as in a localized system g °= exp(— L/^) 11201 . Here £ a o = y/Dtft a is 
the diffusive absorption length, | is the localization length, x a is the ballistic absorption time, 
and g is the total transmission through random medium under diffuse illumination 1251 . When 
the absorption rate 1 jx a exceeds the radiative decay rate through the boundaries of the random 
medium, the resonant tunneling through the modes of the random medium ED is suppressed. 
Because the radiative decay rate diminishes exponentially with an increase of the system size 
L, the amount of absorption needed to suppress the resonant tunneling becomes small. 

In this work, we show that the applicability of SCT in the localized regime is restored when 
sufficiently strong (^„o < L) absorption is added. We also study the size scaling of the position- 
dependent diffusion coefficient and show that in presence of absorption it saturates at a non-zero 
value. This suggests that that the wave transport in absorbing random media can be viewed as 
diffusion, albeit with a renormalized diffusion coefficient, even when localization corrections 
are present and strong. 

2. Model description 

2.1. Motivation for the choice of the model 

The goal of this work is to study the effect of an interplay between absorption and localiza- 
tion on the position-dependent diffusion coefficient. To achieve this goal we selected the model 
geometry based on the following considerations. Because the diffusion coefficient is modified 
due to localization corrections, the system must be in the localization regime in the limit of 
infinitely large system size. Scaling theory of localization predicts [22] that this condition is 
met in systems with physical dimension less than or equal to two. In three dimensional random 
media the disorder strength has to be sufficiently high to reach localization lfT31 . From a com- 
putational standpoint, it is desirable to consider a system with the least dimension. Although a 
one-dimensional (ID) random medium, e.g. a random stack of dielectric slabs l23l or a single 
mode random waveguide l24l . exhibits the phenomenon of localization it is not suitable for 



our study. This is because wave transport in ID shows either ballistic or localized behavior; 
diffusion is not applicable in any parameter range. Instead, we perform numerical simulations 
for 2D disordered multimode waveguides. Because the volume of the system increases as the 
first power of the waveguide length L (for fixed width W), this geometry is often referred to as 
a quasi-lD 12511261. However, unlike ID, random waveguides exhibit both a diffusive behavior 
for £ < L < £, and localization for L > £, . £ denotes the transport mean free path. Furthermore, 
as we argue in Sec. [5] the 2D random waveguides can be fabricated experimentally [27) to test 
the predictions of this work. 

2.2. Position-dependent diffusion coefficient 

To investigate the applicability of a diffusive description of wave transport, we perform the 
ab-initio numerical simulations where all wave interference effects are accounted for without 
any approximations. We consider a scalar, monochromatic wave E(r)e~' m propagating in a 
2D volume-disordered waveguide of width W and length L ^> W that supports N propagating 
modes. The wave field E(r) obeys the 2D Helmholtz equation: 



Here k = co/c is the wavenumber, <5e(r) = (1 + ia)8e r (r), where 8e r (r) is the randomly fluc- 
tuating part of the dielectric constant, and a > is the strength of absorption. The system is 
excited from the left by illuminating the waveguide with N unit fluxes and the wave field E(r) 
throughout the volume of the random medium is computed with a transfer matrix method for 
a given realization of disorder |fT7| . From E(r) we calculate the cross section-averaged energy 
density W(z) and the longitudinal component of flux J z (z). Here z is the axis along the length 
of the waveguide. Local values of the energy density and flux formally define via Fick's law the 
diffusion coefficient D(z) which, in general, may be position-dependent: 



The averages (...) are taken over a statistical ensemble of 10 6 disorder realizations simulated 
on a supercomputer. 

2.3. Dq calculation 

In order to compare our numerical results for D(z) with SCT without fitting parameters, we 
need to obtain the value of the diffusion coefficient unrenormalized by the wave interference 
effects Do = v£/2. Here v is the diffusive speed and £ is the transport mean free path. 

To obtain the value of the transport mean free path in our model we perform a set of simula- 
tions for different waveguide lengths, exploring both the regime of classical diffusion L < £, and 
that of Anderson localization L > <!;, where B, = (n/2)N£. The dependencies of the ensemble- 
averaged conductance (g) and its variance var(g) on the parameter go = (n/2)N£/ (L + 2zo) are 
fitted ifTTl by the analytic expressions obtained by Mirlin in Ref. [28 ] using the supersymmetry 
approach. Here, go the diffusive expression for the average conductance which has the property 
go ~ (g) for L <C t, . The expression for go has two parameters £ and the so-called extrapolation 
length zo EJthat are used to fit numerical data. The presence of scatterers inside the waveg- 
uide causes the effective refractive index in the disordered region to be different from that in 
the empty regions. The refractive index mismatch between can cause surface reflections for the 
waves trying to escape from the system J2]. In the absence of such reflections, zo is related to the 
transport mean free path as zo = {k/A)1. Assuming this relationship between zo and £ we obtain 
a good fit. However, by adjusting I on the order of 2% and treating zo as an independent fitting 
parameter, we achieve even better agreement between the position dependent diffusion coef- 
ficient computed from the ab-initio numerical simulations and from self-consistent theory as 



{V 2 +fc 2 [l + 5e(r)]}£(r)=0. 



(1) 




(2) 



reported in Sec. |4]below. This suggests the existence of some surface reflection in the disorder 
model ifTTl we use in numerical simulations. 

To find the diffusive speed v we use the definition of diffusive flux resolved with respect to 
the direction of propagation [T) 

(Jt\z)) = {v/n)(W{z)) T (D(z)/2)d(W(z))/dz, (3) 

where (z) represent the forward (plus) and backward (minus) propagating components of 
the flux. Combining the two components we find 

2 (W(z)) ' 1 ' 

2.4. % a calculation 

The characteristic absorption time x a can be determined numerically using the condition of flux 
continuity d (J z (z)) jdz = (^(z)) /V The desired diffusive absorption length <i;„o = \/DqT c , can 
be obtained by the proper choice of a in Eq. ([T). 

3. Self-consistent theory 

Self-consistent theory is developed starting with the Green's function G(r,r') of Eq. (HJ 
with 5e(r) = Se r {r) + ia. The disorder-averaged function C(r,r') = (47zWD /cL)(\G(r,r')\ 2 ) 
obeys self-consistent equations in a dimensionless form ll9l [T6l : 



L\ 2 d d 



c(C,C') = 5(C-C'), (5) 



where d(Q = D(Q/Dq and all position-dependent quantities are functions of the longitudinal 
coordinate £ = z/L in the considered waveguide geometry. The quantity C(£, Q which renor- 
malizes the diffusion coefficient is proportional to the return probability at £. Hence, Eq. (O 
formally accounts for wave interference of the time-reversed paths. Eqs. (1516b are solved with 
the following boundary conditions: 

c(C,C')T^(C)^c(C,C') = o (7) 

at £ = and £ = 1. We obtain the numerical solution of Eqs. (|5]|7]i via an iterative procedure 
with an initial guess d(Q = 1. The convergence is usually achieved within five iteration steps. 

4. Analysis of numerical results 

4.1. Interplay between localization and absorption 

First, we consider the effect of localization on the local diffusion coefficient as defined by 
Eq. (O. Fig. QJ compares the ab-initio D{z) (blue curves) and the prediction of SCT from 
Eqs. (|5]|7]i (red curves) in passive systems (no absorption). The behavior of the systems is de- 
termined by the ratio of system length L to localization length E,. The five samples shown 
in Fig. OJ correspond to L/B, — 0.4, 0.7, 1.5, 3.0 and 4.4 covering both L < B, and L > B, 
regimes. We observe that without absorption SCT agrees well with the ab-initio simulations. 
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Fig. 1. Comparison between position-dependent diffusion coefficient d(z) = D{z)/Dq 
found from ab-initio numerical simulations (c.f. Eq. ((2), blue curves) and self-consistent 
theory (Eqs. 1 15 171 , red curves), (a) In passive systems (L/t, = 0.4, 0.7, 1.5, 3.0 and 4.4) 
the diffusion coefficient diminishes in the interior of the system due to the enhanced re- 
turn probability and, hence, stronger localization correction caused by wave interference, 
(b) For fixed length (Z,/| = 4.4), an increase of absorption suppresses the localization 
corrections to the position-dependent diffusion coefficient. The five curves correspond to 
L/^M = 0.0, 3.3, 5.7, 9.8 and 21. (c) When absorption (|/f fl0 = 1.3) is added to the 
five samples shown in (a), the position-dependent diffusion coefficient no longer decreases 
below its saturated plateau value D p , see Eq. (8). In all cases SCT agrees well with the 
ab-initio simulation of wave transport in disordered waveguides. The cause of the small 
deviations are discussed in the text. 



The agreement is achieved for the same values of parameters £ and zo found via the procedure 
described in Sec. [2] We would like to comment on two kinds of systematic discrepancies be- 
tween the numerical data and SCT. (i) We attribute the discrepancy in the boundary regions 
< z < I and L — I <z < Lio the failure of the diffusion approximation at the boundaries. 
This is a well known effect which requires a more sophisticated description, such as that by the 
Milne equation (29], even without localization corrections, (ii) The second discrepancy arises 
at z ~ L/2 deep in the localization regime L^> t,. In the most localized system (L/| = 4.4) 
accessible to us computationally, we find that SCT systematically underestimates the value of 
the position-dependent diffusion coefficient fl6l . An explanation was provided in the work of 
Tian et al lfl8l who pointed out the inability of SCT to describe the effects of resonant tunneling 
which becomes the dominant mechanism of wave transport in a localized system. Based on our 
simulation, we find significant deviations in systems where D(z)/Dq falls below the value of 
~ 0.2, which also agree with the numerical simulations reported in Ref. JTS). 

Absorption has long been considered a nuisance as it tends to both mask and suppress the 
effects of wave localization [20 .26 30, 3 Q. Indeed, absorption suppresses the contribution of 
long propagation paths to the return probability which are responsible for localization effects. 
In Fig.QJ we observe an increase of D(z)/Do — > 1 with a decrease of £, a o (stronger absorp- 
tion) while the sample size L is fixed. The suppression occurs first in the middle of the sample 
(farthermost from the boundaries) where the longest paths are the most probable. We note 
the development of a plateau D(z) ~ D p in the spatial region <i;„o < z < L — ^ ( ,q. This can be 
understood based on the following argument. Due to the presence of absorption, a wave that 
originates at position z in the plateau region can explore only the interval from z — ^ a Q to z + % a o 
which does not include a boundary. Hence, as far as the return probability in Eq. © is con- 
cerned, the effective system size ||26l L e ff ~ 2 x £ a = 2 x y/D p T a , This argument can be made 
more quantitative by neglecting the dd(Q/d£ term in Eq. ©. We solve the self-consistent 
Eqs. (|5]|7|i under the above assumption and obtain the following self-contained equation on 
d p = D p /Do: 

d- l = l + (U/^)d p l/ \ (8) 
which, in the limit of strong absorptions %oial% - > 0, yields an even simpler result d p ~ 
[1 + <?«o/^] _1 - In this expression the parameter g e ff = B, /B, a Q controls the extent of renormal- 
ization of the diffusion coefficient D p , similar to the role played by the conductance go ~ B, /L 
in passive systems Ifl6l . 

A careful comparison between between the numerical results and SCT in Fig. QJ) shows 
a new kind of discrepancy which arises in the localized systems with weak absorption. We 
observe that unlike SCT, the ab-initio D(z) becomes asymmetric with respect to the middle 
of the sample D(z) ^ D(L — z). We hypothesize that this occurs because the absorption tends 
to affect the energy density inside a localized system differently depending on the position of 
the localization center in a transmission measurement [32 ] . The detailed study of these weakly 
absorbing systems goes beyond the scope of this work and will be reported elsewhere. 

4.2. Minimum diffusion coefficient in absorbing systems 

As seen in Fig. QJ, in the passive systems the value of the diffusion coefficient in the middle 
of the sample D(z = L/2) approaches zero with an increase of L. As the boundaries are further 
removed from z = L/2, it becomes progressively more difficult for a wave to escape and the 
return probability keeps increasing. As discussed above, the same argument no longer applies in 
an absorbing random medium. In Fig. QJ we plot the position-dependent diffusion coefficient 
in absorbing (£,/£, a o = 1.3) systems of different size. Unlike passive systems in Fig. QJ, the 
diffusion coefficient stops decreasing after the system size L exceeds the length ~ 2£, a o and 
develops a plateau where it is no longer sensitive to the boundaries of the system. 
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Fig. 2. Existence of the minimum diffusion coefficient is seen from the evolution of 
d^2 = D(z = L/2)/Dq with the increase of the system size. In the passive systems 
(if; /£, a Q = 0, open circles) the limit is expected to be zero. In absorbing systems (£, /^,o = 
0.3, 0.7, 1.3, 2.9 shown as cross, diamond, upward and downward triangle symbols re- 
spectively) saturation corresponds to formation of the plateau region seen in Figs. [Tp,c. 
The saturation value D p increases monotonically with an increase of £,/t, a o- Five solid 
lines are obtained from the self-consistent theory Eqs. (1517b for each value of the absorp- 
tion strength. Qualitative prediction of the minimum value of position-dependent diffusion 
coefficient in SCT is supported by the numerical simulations. The agreement is also quan- 
titative for D p /D > 0.2. 

Fig. [2]shows saturation of the diffusion coefficient in the middle of the sample D(z = L/2) 
with an increase of the sample size L. Importantly, the saturated value D p is determined by the 
ratio £, /£, a Q, see Eq. ©. As we discussed above, SCT underestimates the value of the diffusion 
coefficient when the latter becomes small. However, the phenomenon of saturation is observed 
in both ab-initio simulations and SCT. 

4.3. Universal scaling of position-dependent diffusion 

Single parameter scaling 11221 is a key concept in mesoscopic physics. It predicts that when the 
size of the system increases, the evolution of the average conductance, its variance 11331 as well 
as its entire distribution ll34l are determined by the average conductance. In other words, the 
mesoscopic electronic transport through a disordered system is universal l35l - independent 
of the microscopic parameters of disorder. The concept of universality can be directly applied 
also to the electromagnetic wave transport 0. Recently, it was shown that in the absence of 
absorption the concept of position-dependent diffusion coefficient is also universal |[T6HT8l - 
the shape of D(z) is uniquely determined by the value of conductance. 

In absorbing systems, there appears the second parameter related to the absorption length 
scale [23 , 26 36 37 1 . Throughout this work we use £, /£, a o as such a parameter. The agreement 
between the numerical simulations and SCT demonstrates that the position-dependent diffu- 
sion coefficient is uniquely determined by two dimensionless parameters go ~ £, /L (related to 
conductance in passive systems) and ^/^ a Q - the only parameters in self-consistent Eqs. (07]i. 
The parameter zo entering into Eq. (O is not essential as it affects the behavior of D(z) only in 
the boundary regions where its behavior is not universal |fP7l . 

Eq. ^ describing the plateau value of the diffusion coefficient in absorbing systems depends 
only on one parameter ^ /^„o- This suggests a new universal behavior. In Fig. [3] we compare 





1 


1 

=2.9 


■ \r~~ — v v — 


v — 




=1.3 


\\>~~"~" — — a 


0^ 


=0.7 




_ *au 


=0.3 
=0.0 





+ - L/t =0.4 N=20 

^> -L/£=0.6N=20 

tf - L/i =0.8 N=20 

* - L/C =1.1 N=20 



11 0.1 1 10 

f/£o0 X [l - d 1/2 (oo)] 

Fig. 3. Universal scaling of the local diffusion coefficient in the middle of the sample = 
D(z = L/2)/Dq is described by a single parameter. The ab-initio numerical simulations 
in disordered waveguides with absorption (symbols) agree well with Eq. {§}, which we 
derived based on SCT. 

this prediction of SCT (solid line) to our numerical results in different disordered absorbing 
waveguides (symbols) and observe a good agreement. 

5. Conclusions 

In this work we compared the results of the ab-initio numerical simulation of wave trans- 
port through disordered absorbing waveguides to the prediction of self-consistent theory with 
position-dependent diffusion. This theory has the power to provide a simple intuitive descrip- 
tion of an intricate interplay between localization effects due to wave interference and absorp- 
tion. The following results were obtained, (i) We demonstrated that in localized regime (L > t, ) 
for weak/strong absorption SCT provides the correct qualitative/quantitative description of the 
renormalizations in position-dependent diffusion coefficient, (ii) We pointed out the key dis- 
tinction in scaling of the position-dependent diffusion in passive and absorbing media. In the 
former, the diffusion coefficient vanishes in the limit of infinite size. In contrast, in absorbing 
random media, the local diffusion coefficient always remains finite so that true localization in 
the sense of vanishing diffusion is never reached. This conclusion appear to be borne out also 
by supersymmetric field theory [38 1 ■ 

We would like to point out that a similar (but not equivalent) effect occurs in the mesoscopic 
electronic systems due to the effects of dephasing 1291 which introduces the cutoff length scale 
beyond which waves no longer contribute coherently to the return probability. Instead, these in- 
elastically scattered (de Broglie) waves of the electrons contribute to the incoherent background 
which only obscures the localization phenomena. In contrast, optical absorption removes pho- 
tons while the remaining component retains all information about interferences. This observa- 
tion points to the possibility of observation of the phenomenon of position-dependent diffusion 
in optical system that we discuss next. 

There are three main conditions that would allow the direct experimental observation 
position-dependent diffusion. First, without absorption the system must be in the localization 
regime in the limit of infinitely large size. Second, the absorption, or more generally loss, can- 
not be too strong so that the localization effects are not completely washed out. Third, because 
the renormalization of diffusion occurs away from the boundaries, the interior of the medium 
has to be experimentally (non-invasively) monitored. These conditions can be met in the two 
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dimensional disordered optical waveguides where the position-dependent diffusion is indeed 
observed 11271 . 
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